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It is demonstrated that using the Fridkin and Gardon - Narayanaswami algorithms makes it possible to de¬ 
scribe conditions of glass hardening with a sufficient accuracy and provides enough reason to calculate non¬ 
linear glass treatment regimes, which makes it possible to predict new consumer properties in sheet glass. 


One of the most important operations in thermal treat¬ 
ment of sheet glass is hardening. Articles treated in this way 
have high strength parameters and a safe type of destruction. 
However, they cannot undergo machine treatment (cutting, 
drilling, grinding, etc.) due to their self-sustaining destruc¬ 
tion. It is essential to develop such method for thermal treat¬ 
ment of glass that would yield a product combining the pos¬ 
sibility of mechanical treatment with a high level of strength. 

The only solution so far found for this problem is com¬ 
plex thermal treatment (CTT) of sheet glass [1], This 
method, similarly to traditional hardening, can be implemen¬ 
tation both on horizontal and on vertical hardening lines with 
a minimum correction required of the technological process. 
The only problem encountered by production engineers is 
the absence of algorithms for estimating the properties of 
sheet glass after CTT. 

The method for calculating the process of hardening un¬ 
der a constant intensity of heat exchange with the coolant 
(a = const) is sufficiently well investigated [2]. 

At the same time, implementation of nonlinear heat treat¬ 
ment conditions with a variable heat transfer coefficient in 
cooling (a ^ const) would make it possible to significantly 
expand the range of consumer properties of sheet glass, 
which would open additional prospects in producing lami¬ 
nated protective glass with a significantly decreased weight 
of articles while preserving the safety parameters. Further¬ 
more, it has been established in preliminary studies that the 
CTT process is going to be less energy-consuming than clas¬ 
sical hardening. 

The proposed algorithm is based on the studies [1] mak¬ 
ing it possible to calculate temporary and residual stresses 
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taking into account the nonlinear heat transfer coefficient a 
throughout the entire thermal treatment. 

The difficulties in the development of a method for cal¬ 
culating stresses in CTT technology are related to the impos¬ 
sibility of using the classical Fourier solutions for the heat 
transfer coefficient that is nonlinear in time. Using the Fou¬ 
rier method [3] in these conditions, a researcher is confronted 
with objective inadequacy of the obtained solutions. 

Let us consider cooling of glass for a time period 
ij = 20 sec with an air flow a, and subsequent exposure of 
glass under free convection in the context of a qualitative es¬ 
timate (Fig. 1). During intense cooling the temperatures of 
the outer and inner layers (curves 1 and 2 ) correlate with the 
true temperature field. However, if the Fourier equations are 
used when the cooling regime changes (points A and B ), the 
temperature of the outer layer “freezes” (curve 3 ) and the 
outer layer temperature abruptly drops (curve 4 ), which con¬ 
tradicts common sense. 

In this case balance equations are preferable. An exam¬ 
ple of such equations was used by R. Z. Fridkin and 
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Fig. 1 . Temperature fields in the course of complex thennal treat¬ 
ment (calculated by the Fourier method). 
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Fig. 2. Temperature fields in the course of complex thermal treat¬ 
ment (calculated by the Fridkin method). 


The glass surface temperature can be calculated in two 
stages: 

in the first approximation 

T' s (x 0 , x + At) = 1.875r(x 1 , t) + 

0.375r(x 3 , t)- 1.25r(x 2 , x), 

and finally, taking into account surface effects: 

7’ s (x 0 ,t + At) = 

ATjx^ + aJ.-z^EiTD-O.lEiT^)] 
a a + A 


O. V. Mazurin [4] in calculating heating of glass. That algo¬ 
rithm implies the calculation not only of the conductive com¬ 
ponent of heat transfer but of the radiant component as well, 
since at low cooling intensity the latter constitutes up to 30% 
[5]. The distinctive feature of the algorithm described in the 
present paper is calculation of the integral radiation compo¬ 
nent using the Stefan - Boltzmann formula (without using 
power series), which makes it less cumbersome. 

The qualitative representation (Fig. 2) of variation of the 
temperatures of the inner (curves 1 and 3 ) and outer (curves 
2 and 4 ) layers yields a more sensible result in the form of a 
temporary “outburst” of the outer layer temperature (curve 
2 ) due to the conductive component (curve 4 ) from the inner 
layers, and the calculation algorithm is as follows. 

The temperature of the /th layer of glass at time x + Ax is 
equal to 

T (x,-, x + Ax) = T (x,., x) + — [<? r (x,-, t) + q c (x ; , x)], 
cp 


where a a is the convective heat transfer; s s is the reduced de¬ 
gree of blackness of the cooling surfaces; A = 2X*/{Q.\d ). 

The components E(T' % ) and E ( T) both represent inte¬ 
gral radiation of the article at the surface temperature and at 
the ambient air temperature, respectively. They are deter¬ 
mined from the Stefan - Boltzmann formula [6]: 


E(T) = 0.37413 x 10 9 £- 


AX 


*=°X 5 


exp 


14,410 A 
XT J 


-1 


where X is the wavelength. 

Further calculations are perfonned on the basis of the 
points with coordinates 


d 
2 n 


(n-i) 


where T(x [ , x) is the temperature at the preceding moment; 
c is the specific heat; p is the glass density; q x and q c are the 
radiant and conduction components of specific heat at the 
preceding time moment. 

The sheet thickness d is divided into 2 n layers. Subjected 
to the condition of a symmetrical temperature distribution 
with respect to the central section of the plate, the coordinate 
of the /th layer (except for the surface layer) is equal to 

X,- = -p- (2/ - 1), 

4 n 

and it is sufficient to set n = 5. 

The conductive component of any layer except for the 
first one amounts to 


q c (x,., x) = X* 


r(x,- 1 ,x)-2r(x,-,x) + r(x, + i,x) 

(0.W) 2 


and for the first layer 


with simultaneous conversion to temperature values in °C: 
for layers (/ = 1 ... n - 1) 


/(x',x) : 


T ( X n -b'C j )+T(X n - i+ 1,T;) 


-273.15; 


for the boundary layers (z = 0, n ) 

t(x\,Xj ) = t{x n _ i ,x j )- 273.15. 

Next, the real viscosity r) of glass is determined using the 
Tamman - Fulcher and the Arrhenius formulas: 


T) (x,., x ) = exp 


B 


v T{x i ,x j )-C J 


at 9 < log q < 13; 


r| (x ; ., %j ) = 10 T( ' x ‘ ,T/ * at 9 > log r| > 13, 


q c (X[, x) = X* 


2T s (x 0 ,x)-3T(x 1 ,t)+ T(x 2 ,t) 
(0 Id) 2 


where X* is the effective thermal conductivity of glass. 


where A, B, C, A\ and B' are constants depending on the 
glass composition. 

Subsequent calculation of temporary and residual stresses 
corresponds to the procedure described in [1], 
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Fig. 3. Temporary stresses in sheet glass during hardening. 

Based on the developed complex algorithm, a program 
was developed and a reference calculation was performed 
based on experimental data of R. Gardon [5] with the initial 
conditions; t(0) = 648°C, a = 222.6 W/(m 2 • K), and 
d= 6.1 mm. The correlation results are shown in Fig. 3. The 
dots indicate experimental data and the lines indicate tempo¬ 
rary stresses calculated using the proposed algorithm. 

The mean deviation of estimated results from the experi¬ 
mental data amounted to 9.68%. 

Based on the considered method, CTT regimes were de¬ 
veloped and sheet glass was hardened in accordance with 


these regimes. The resulting glass has high mechanical 
strength and is susceptible to mechanical treatment. 

Thus, using the Fridkin and Gardon - Narayanaswami 
algorithms makes it possible to describe conditions of glass 
hardening with a sufficient accuracy and provides enough 
reason to calculate nonlinear glass treatment regimes, which 
makes it possible to predict new consumer properties in sheet 
glass. 
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